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Abstract 


We report on the observations of the quasar NRAO 530 with the Event Horizon Telescope (EHT) on 2017 April 5—7, 
when NRAO 530 was used as a calibrator for the EHT observations of Sagittarius A". At z = 0.902, this is the most 
distant object imaged by the EHT so far. We reconstruct the first images of the source at 230 GHz, at an unprecedented 
angular resolution of ~20 pas, both in total intensity and in linear polarization (LP). We do not detect source variability, 
allowing us to represent the whole data set with static images. The images reveal a bright feature located on the southern 
end of the jet, which we associate with the core. The feature is linearly polarized, with a fractional polarization of ~5%-— 
8%, and it has a substructure consisting of two components. Their observed brightness temperature suggests that the 
energy density of the jet is dominated by the magnetic field. The jet extends over 60 pas along a position angle ~ —28°. 
It includes two features with orthogonal directions of polarization (electric vector position angle), parallel and 
perpendicular to the jet axis, consistent with a helical structure of the magnetic field in the jet. The outermost feature has 
a particularly high degree of LP, suggestive of a nearly uniform magnetic field. Future EHT observations will probe the 
variability of the jet structure on microarcsecond scales, while simultaneous multiwavelength monitoring will provide 


insight into the high-energy emission origin. 


Unified Astronomy Thesaurus concepts: Extragalactic astronomy (506) 


1. Introduction 


NRAO 530 (1730—130, J1733—1304) is a flat-spectrum 
radio quasar (FSRQ) that belongs to the class of bright y-ray 
blazars and shows significant variability across the entire 
electromagnetic spectrum. The source was monitored by the 
University of Michigan Radio Observatory at 4.8, 8.4, and 
14.5 GHz for several decades'*! until 2012. The quasar 
underwent a dramatic radio outburst in 1997 (Aller et al. 
2009), during which its flux density at 14.5GHz exceeded 
10Jy, while the average value is ~2Jy. Since 2002, 
NRAO 530 has been monitored by the Submillimeter Array 
(SMA; Maunakea, Hawaii) at 1.3 mm and 870 um. 15? The data 
show high amplitude variability at 1.3 mm from 1—4 Jy, with 
the brightest outburst in the end of 2010 contemporaneous 
with prominent y-ray activity (Williamson et al. 2014). 
Ramakrishnan et al. (2015) found a statistically significant 
correlation between the radio and y-ray light curves of 
NRAO 530, based on measurements at 37GHz by the 
Metsdhovi Radio Observatory (Aalto University, Finland) and 
y-ray fluxes measured with the Large Area Telescope (LAT) on 
board the Fermi Space Gamma-Ray Telescope. The blazar is 
included in a sample of active galactic nuclei (AGNs) observed 
with the POLAMI program (Agudo et al. 2014), which 
provides both the flux density and polarization measurements 
at 3mm and 1.3 mm started in 2009. According to the data 
during the perion 2007-2021, the degree of linear polarization 
(LP) at 1.3mm changes from 1% to 15%, with circular 
polarization detections at a level of 196—296 (private 
communication). 

The quasar is known as a violently variable object at optical 
wavelengths. According to the SMARTS optical/IR monitoring 
program of blazars (Bonning et al. 2012), NRAO 530 has a 
fractional variability amplitude of ~68% and ~35% at optical and 
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IR wavelengths, respectively, with an average magnitude of ~14 
at the K band (\gp~2.2m) and ~17.5 in the R 
band (Acpp~ 658 nm). Foschini et al. (2006) reported a 
short hard X-ray flare detected serendipitously by the IBIS/ 
ISGRI detector on board INTEGRAL on 2004 February 17, 
although the timescale of the flare, «1hr, questions the 
association with the FSRQ. However, a moderate increase of 
the LP at 2 cm contemporaneous with the flare was observed as 
well. The recent 10 yr Fermi LAT catalog, 4FGL (Abdollahi et al. 
2020), lists the average y-ray flux of the source at 0.1—300 GeV as 
(1.8 + 1.2) x 1077 photcm ? s™". During this period NRAO 530 
underwent the highest amplitude ^-ray outburst in 2010 October- 
November, when the y-ray flux rose to 10 ^ photcm ?s^! 
(Williamson et al. 2014). This y-ray activity coincided with a 
significant brightening at optical wavelengths and an increase of 
the optical linear degree of polarization up to 1506. >? 

NRAO 530 possesses a highly relativistic jet. The source is 
intensively monitored with the Very Long Baseline Array 
(VLBA) at 15 GHz within the MOJAVE program (e.g., Lister 
et al. 2016) and at 43 GHz within the VLBA-BU-BLAZAR 
program (Jorstad et al. 2017; Weaver et al. 2022). VLBA 
images show a radio jet dominated by the core located at the 
southern end of the jet. At 15 GHz, the jet extends up to 10 mas 
to the north, with a weak feature observed ~25 mas from the 
core. Based on 11 moving features, Lister et al. (2019) reported 
a median apparent speed of 12.3 +0.6c, with a maximum 
apparent speed of 27.3 + 1.0c. At 43 GHz, the jet contains a 
stationary feature located ~0.3 mas from the core and curves to 
the northwest —0.8 mas from the core. Eight moving knots 
detected at 43 GHz over 2007-2018 exhibit a wide range of 
apparent speeds, from 2.6 + 0.5 c to 37.2 + 0.2 c (Weaver et al. 
2022), with three knots revealing a strong acceleration after 
passing the stationary feature, while knot B3 ejected in the 
beginning of 2010 demonstrates both accelerations from 4.1 c 
to 7.7 c and to 33.7 c and then a deceleration to 8.5 c. Using 
properties of the light curves of moving knots and their 
apparent speeds, Jorstad et al. (2017) estimated the average 
Lorentz and Doppler factors of the jet to be I — 8 and ~ 9, 
respectively, with a jet viewing angle of ~3°. These values of 
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the Lorentz and Doppler factors are close to the minimum 
values obtained by Liodakis et al. (2018) for the jet of 
NRAO 530. These authors estimated the variability of the 
Doppler factor by analyzing the 15 GHz light curve provided 
by the Owens Valley 40m Radio Telescope after taking into 
consideration brightness temperature arguments and the 
apparent speeds of knots at 15GHz; this resulted in 
Tmax ~ 40 and max ~ 21. 

On kiloparsec scales, the quasar shows a weak two-sided jet 
elongated in the east-west direction over 2”5, with a very 
diffuse eastern part and a bright, compact knot on the western 
side, located about 1” from the core (Kharb et al. 2010). The 
kiloparsec-scale jet is almost perpendicular to the parsec-scale 
jet and has a complex radio morphology: the eastern part is 
similar to that expected for Fanaroff—Riley I (FRI) type radio 
galaxies, while the western part can be classified as FRII type. 

NRAO 530 has a redshift of z — 0.902 (Junkkarinen 1984), 
for which 100 pas corresponds to a linear distance of 0.803 pc 
(Planck Collaboration et al. 2016, Hyp = 67.7 km s ! Mpc |, 
Qm = 0.307, and €, = 0.693). The source contains a super- 
massive black hole, the mass of which is currently uncertain, 
with estimates ranging from 3x 105 M (Keck 2019) to 
2 x 10? M, (Liang & Liu 2003). 

Due to its brightness, compactness, and close location on the 
sky to the Galactic Center, NRAO 530 is frequently used as a 
calibrator in very long baseline interferometry (VLBI) 
observations of the radio source Sagittarius A’ (Sgr A* ) located 
in the center of the Milky Way (e.g., Lu et al. 2011). In this 
paper we report on the VLBI observations of NRAO 530 in 
2017 April, when the source was employed as a calibrator 
during the Event Horizon Telescope (EHT) observation of 
Sgr A* at 1.3mm (EHTC et al. 2022a). During the EHT 
observation, NRAO 530 was in a moderate activity state at J- 
ray energies with the flux ~1 x 107 phot cm ?s ! as well as 
at optical (R-band magnitude w17.8) ”“ and radio wavelengths 
(the flux density at 37 GHz ~3.7Jy).'* According to the 
VLBA monitorings at 43 GHz (Weaver et al. 2022) and 
15 GHz (Lister et al. 2019), no superluminal components were 
ejected in 2017. The EHT observation of NRAO 530 provides 
unique insight into the subparsec scale structure of the quasar 
with microarcsecond resolution, which are important for 
understanding the blazar physics in general and properties of 
the NRAO 530 jet in particular. The outline of the paper is as 
follows: in Section 2 we describe the EHT observations of 
NRAO 530, their reduction and general data set properties, in 
Sections 3 and 4 we present results of total intensity and LP 
imaging, respectively, in Section 5 we discuss our results, and 
in Section 6 we summarize our findings. 


2. Observations and Data Calibration 


The EHT observed NRAO 530 on three consecutive nights 
on April 5-7 during the 2017 campaign (MJD 57848-57850), 
as a calibrator of Sgr A (EHTC et al. 2022a, 2022b). The 
observations were performed with the full EHT 2017 array of 
eight telescopes located at six geographical sites: the Atacama 
Large Millimeter/submillimeter Array (ALMA) and the 
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Atacama Pathfinder Experiment (APEX) telescope in Chile; the 
Large Millimeter Telescope Alfonso Serrano (LMT) in 
Mexico; the IRAM 30 m telescope (PV, joined only on April 
7) in Spain; the Submillimeter Telescope (SMT) in Arizona; 
the James Clerk Maxwell Telescope (JCMT) and SMA in 
Hawaii; and the South Pole Telescope (SPT) in Antarctica. 
Two 2GHz wide frequency bands, centered at 227.1 GHz 
(LO), and 229.1 GHz (H1), were recorded. The observations 
were carried out using dual feeds, right-hand and left-hand 
circularly polarized (RCP ad LCP, respectively), for all stations 
other than ALMA and JCMT. ALMA recorded dual LP, which 
was subsequently converted at the correlation stage to a circular 
basis by PolConvert (Marti-Vidal et al. 2016; Goddi et al. 
2019). The JCMT only observed the RCP component, which 
we utilized to approximate total intensity (this approximation is 
correct as long as the fractional circular polarization of the 
source can be neglected), and we typically omit JCMT 
baselines for the polarimetric imaging. The EHT array setup 
is detailed in EHTC et al. (20192, 2019b). 

Recorded signals were correlated at the MIT Haystack 
Observatory and the Max-Planck-Institut für Radioastronomie, 
Bonn. Subsequent data reduction procedures are described in 
EHTC et al. (20190), Blackburn et al. (2019), and Janssen et al. 
(2019). There were minor updates to the calibration with 
respect to the EHT results published earlier (EHTC et al. 
2019a; Kim et al. 2020), particularly regarding the telescope 
sensitivity estimates and complex polarimetric gains calibra- 
tion. These updates are identical to those described in the EHT 
Sgr A* publications (EHTC et al. 2022a, 2022b). Flux density 
on the short intra-site baselines (ALMA-APEX and SMA- 
JCMT) was gain-calibrated to the simultaneous ALMA-only 
flux density of 1.6Jy, reported by Goddi et al. (2021), 
providing an additional constraint on the amplitude gains for 
the stations with a colocated partner (network calibration; 
Blackburn et al. 2019). The polarimetric leakage calibration 
follows procedures outlined in EHTC et al. (2021), where 
D-terms of stations with an intra-site VLBI baseline (ALMA, 
APEX, SMA, and JCMT) were calculated through a multi- 
source fitting procedure (Martí-Vidal et al. 2021), and several 
estimates for D-terms of remaining stations were reported. Here 
we employ the fiducial D-terms given by Issaoun et al. (2022) 
based on the analysis presented in EHTC et al. (2021). Some 
imaging algorithms that we use ignore the leakage calibration 
and estimate D-terms independently, providing an additional 
consistency check of the calibration procedures and their 
impact on the resulting images (see Section 3). 

The (u, v)-coverage of the EHT observations of NRAO 530 
is given in Figure 1. We show the coverage on each individual 
day, as well as aggregated over the entire observing campaign. 
The longest projected baseline reaches 8.66 GA for SPT-PV, 
which determines the instrumental angular resolution of 24 pas, 
measured as the minimum fringe spacing. Inspection of the 
data revealed no indication of the source structure evolution 
during the observations spanning only about 52 hr on 2017 
April 5—7, with the total on-source integration time of 76 
minutes (see Figures 2—3). For that reason, the entire data set 
was used simultaneously for static imaging and model fitting. 

The NRAO 530 data set self-consistency has been rigorously 
verified as a validation of the EHT Sgr A* results, presented in 
EHTC et al. (2022b). The signal-to-noise ratio (S/N) of 
individual NRAO 530 detections (band-averaged in frequency, 
scan-averaged in time) typically exceeds 100 on baselines to 
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Figure 1. EHT (u, v)-coverage of the NRAO 530 observations on 2017 April 5, 6, and 7, and all days aggregated. Each colored point corresponds to a single VLBI 
scan of 3-4 minutes. ALMA participated in the observations on 2017 April 6 and 7. Dashed circles indicate the fringe spacing of 50 and 25 pas. “Chile” represents the 


stations ALMA and APEX. “Hawaii” represents the stations SMA and JCMT. 


ALMA, and is about an order of magnitude lower on remaining 
baselines; see the top panel of Figure 2. We estimated an 
additional 2% of complex visibility uncertainties due to variety 
of nonclosing effects (e.g., residual polarimetric leakage, 
frequency, and time smearing); see EHTC et al. (2019c, 
2022b). Inspecting the calibrated visibility amplitudes in the 
middle panel of Figure 2, we notice that the shortest inter-site 
baseline (LMT-SMT) provides a lower limit on the compact 
emission flux density corresponding to about 0.6 Jy, which 
implies the presence of up to | Jy of milliarcsecond-scale flux 
that is resolved out by the EHT but measured on the intra-site 
VLBI baselines and by the ALMA-only observations. Addi- 
tionally, in the bottom panel of Figure 2 we present the 
measured visibility amplitudes as a function of the baseline 
length, projected at the jet PA (see Section 3). We notice that 
the amplitudes on baselines shorter than about 2.5 GA can be 
approximated by a Gaussian of FWHM 55 yas and 0.6 Jy 
amplitude, roughly informing us about the size and the total 
flux density of the observed compact feature. Measurements on 
baselines longer than 2.5 GA indicate a deviation from the 
Gaussian model, in agreement with the resolved source 
consisting of more than one component Both EHT 
Sgr A* calibrators, NRAO 530 and J1924-2914 (Issaoun et al. 
2022) are fairly compact and bright sources at 230 GHz. 
During the week-long EHT observing run, their source 
structure and flux density did not show significant variations 
(see Figure 3). These properties allowed us to use both sources 
for quantifying the effects of the residual errors in the 
characterization of the antenna gains. Mitigating such effects 
was crucial for discerning between the rapid intrinsic variability 
of Sgr A* and the apparent flux density variations produced by 
imperfect calibration (EHTC et al. 2022c, 2022d; Wielgus et al. 
2022). Details about the procedures and methods used to obtain 
the antenna gains from the two calibrators and the subsequent 
transfer to the Sgr A* data sets are reported in Section 3.2 in 
EHTC et al. (2022b). 


3. Total Intensity Imaging 


For the imaging of NRAO S530, we have utilized data 
calibrated by the EHT-HOPS pipeline (Blackburn et al. 2019), 
combined over all days (2017 April 5—7) and both frequency 


bands (227.1 and 229.1 GHz). Nonzero closure phases shown 
in Figure 3 reveal resolved and nontrivial structure of the 
source. We used three analysis methods across five different 
implementations. These are: imaging via inverse modeling (the 
traditional CLEAN method), forward regularized maximum 
likelihood (RML) modeling (eht-imaging and SMILI 
algorithms), which are described in detail in EHTC et al. 
(20194), as well as two posterior exploration methods based on 
a Markov Chain Monte Carlo (MCMC) scheme—a D-term 
Modeling Code (DMC; Pesce 2021) and a dedicated Bayesian 
image reconstruction framework Themis (Broderick et al. 
2020a). For each of the methods, the freedom to decide 
resolution, field of view, or amount of data coherent averaging 
was left to the expert imaging subteams, who made such 
decisions based on the specific properties of each algorithm. As 
an example, MCMC methods tend to average data more to 
reduce the computational complexity of the likelihood evalua- 
tion. Below we give a short description of each method, with 
some details on the particular application to the NRAO 530 
imaging. 


3.1. DIFMAP 


For CLEAN imaging, we have used the DIFMAP software 
package (Shepherd 1997). We have constructed multiple 
manual images prior to the final process to obtain a set of 
CLEAN windows corresponding to the source structure and 
agreeing best with the data. This window set includes a large 
circular Gaussian with a radius of 1 mas, located in the center 
of the map, which models the large-scale extended emission 
present in the source. The standard process of multiple 
consecutive CLEAN iterations, along with self-calibration 
procedures, was employed to achieve the best normalized x? 
values between the data and the model for the closure phases, 
closure amplitudes, and visibilities. We have used data 
averaged in 10s segments, a field of view of 1024 x 1024 
pixels with a pixel size of 2 uas, and uniform weighting to 
optimize the image resolution. 


3.2. eht-imaging 


The eht-imaging code (Chael et al. 2016, 2018) 
implements the RML-based approach to imaging. In this 
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Figure 2. Top: signal-to-noise ratio (S/N) of NRAO 530 detections. Colors 
follow the convention of Figure 1. The entire data set (2017 April 5-7, both 
frequency bands) is shown. The data were averaged coherently in 240s 
intervals. The circles denote primary baselines, while the diamonds denote 
redundant baselines (i.e., obtained with APEX or JCMT, rather than with 
ALMA or SMA). The cluster of S/N > 100 circles corresponds to ALMA 
baselines. Middle: the same as the top panel, but for the correlated flux 
densities (visibility amplitudes) in the a priori calibrated data set. Bottom: 
visibility amplitudes as a function of baseline length projected at the jet 
position angle (PA) of —25 deg (see Section 3). The dashed line represents a 
Gaussian component of FWHM 55 pas and 0.6 Jy amplitude. 


method, an updated image is being proposed at every iteration, 
while the algorithm optimizes a cost function comprising the 
error term describing consistency with the data (predominantly 
through closure phases and log closure amplitudes; Blackburn 
et al. 2020), and regularization terms, similar to Bayesian 
priors. The latter may promote image properties such as 
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smoothness, sparsity, compactness, or consistency with the 
imaging prior—here chosen to be a 60 pas FWHM circular 
Gaussian. For the NRAO 530 analysis, we assumed a field of 
view of 200 uas and a 64 x 64 pixel grid, and used scan- 
averaged data. The algorithm proceeds with an iterative 
sequence of gain calibration and cost function optimization. 


3.3. SMILI 


SMILI (Akiyama et al. 2017) is another RML imaging 
library that reconstructs interferometric images utilizing a 
similar set of regularizers including weighted-/, (£P), total 
variation, total squared variation, and the maximum entropy 
regularizer. Prior to the imaging, we rescale the intra-site 
baseline flux density of NRAO 530 to 0.9Jy to remove the 
contributions from the large-scale extended emission. During 
the imaging, we adopt a field of view of 300 uas and a 
150 x 150 pixel grid, initializing the algorithm with a circular 
Gaussian prior, and using 10 s averaged data set. Subsequently, 
the imaging algorithm minimizes the regularized cost function. 
After 1000 iterations, self-calibration of the stations gains is 
performed, and the image prior is updated. This imaging and 
self-calibration loop is repeated until the reconstruction 
converges. 


3.4. DMC 


DMC formulates the imaging problem in terms of posterior 
exploration, which is carried out using a Hamiltonian Monte 
Carlo sampler as implemented in the PyMC3 Python package 
(Salvatier et al. 2016). The output of DMC is a collection of 
samples from the joint posterior over the image structure and 
calibration quantities. Rather than carrying out an iterative self- 
calibration procedure, DMC fits the complex gains at every 
station simultaneously with the full Stokes image structure. For 
NRAO 530 imaging, a field of view of 300 uas and a 30 x 30 
pixel grid were used, the data were scan-averaged, and D-terms 
were fixed following Issaoun et al. (2022). A comprehensive 
description of the DMC model and implementation can be found 
in Pesce (2021). 


3.5. Themis 


Themis is a framework for fitting VLBI data with 
an MCMC posterior space exploration method (Broderick 
et al. 2020a). For NRAO 530 data analysis, the image was 
modeled with a spline raster defined using 6 x 5 control points, 
and an adaptively selected field of view (Broderick et al. 
2020b). The data were scan-averaged. We solved for D-terms 
with Themis in order to provide a consistency test for the 
fiducial D-terms of EHTC et al. (2021) and Issaoun et al. 
(2022); see Section 4.1. 


3.6. Fiducial Images of NRAO 530 


Figure 4 presents final (fiducial) total intensity images 
obtained by the methods discussed in Sections 3.2—3.5. For the 
presentation purposes, all images are shown using the same 
field of view of 128 x 128 uas, more narrow than the field of 
view typically used for the imaging. The contribution of 
emission outside of this region is relatively small and consistent 
with the image noise level, which, in the case of the EHT, is 
rather high and dominated by the systematic uncertainties 
related to sparse (u, v)-coverage. There is a high degree of 
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Figure 3. NRAO 530 closure phases measured on a selected nontrivial triangles. LO band data points averaged in 240 s intervals are shown. There is no statistically 
significant indication of structural variation between the three consecutive observing days. Predictions of the models obtained with a variety of imaging algorithms, 
discussed in Section 3, are indicated with thin black lines. The blue-shaded region indicates the range of predictions from different reconstructions presented in this 


paper. 


agreement between the different imaging methods in terms of 
the source structure and the relative brightness of individual 
features. All images show the jet to be elongated from the 
southeast to the northwest over ~60 uas. A visual inspection 
suggests that the jet consists of two dominant features (see 
Figure 5), with the brightest one, CO, located at the southeast 
end of the jet and itself exhibiting substructures COa and COb. 
However, there is also a difference between the models— 
DIFMAP and DMC images contain a very low-brightness 
extended feature that turns to the northeast, almost by ~90° to 
the main direction of the jet, the eht- imaging image shows 
a hint of such a feature, while the SMILI and Themis 
methods do not reveal its presence. This discrepancy is likely 
related to the latter two algorithms promoting sparsity and 
compactness in the reconstructed images. We have designated 
this ambiguous feature as C2 in Figure 5, marked by a dashed 
line ellipse. Although the knot appears in the fiducial method- 
averaged image if we increase the dynamic range, its brightness 
peak is on the level of the image noise. In addition, DI FMAP, 
eht-imaging, and SMILI images correspond to very 
similar peak brightness values, but the Themis and DMC 
images indicate somewhat brighter and weaker brightness 
peaks, respectively. 

Figure 3 shows consistency between closure phase measure- 
ments on selected triangles and corresponding predictions of 
the fitted models. Table 1 provides goodness of fit normalized 
x statistics computed for the models shown in Figure 4 for 
self-calibrated visibility amplitudes OG ues closure phases 
Coa and log closure amplitudes hec A). The estimated total 
flux density of the compact source model is also provided. This 
quantity is particularly challenging to constrain with the EHT 
array, which in 2017 lacked baselines between 5 MA and 
500MA projected length, sensitive to milliarcsecond-scale 
emission. Based on the presented statistics, there is no strong 
preference for any particular model. Therefore, we have 
combined the final images constructed by the different 
methods, shown in Figure 4, to produce afiducial method- 
averaged total intensity image of NRAO 530, presented in 
Figure 5. The averaging procedure involves aligning images to 
maximize their cross-correlation, re-gridding them to a 
common resolution and field of view, and taking a standard 
arithmetic mean on a pixel-by-pixel basis. It follows the 


practice of EHTC et al. (2019d, 2022c), Kim et al. (2020), and 
Issaoun et al. (2022), and it is designed to make the resulting 
image-domain morphology more robust against method- 
specific systematics. 


3.7. Total Intensity Component Modeling 


We have modeled the visibilities, calibrated to the fiducial 
image, with circular Gaussian components using the task 
modelfit in DIFMAP. We also fitted a Gaussian model to the 
data in the eht - imaging framework, using closure-only data 
products to constrain the resolved source structure. In the latter 
case, the simplest model representing data properties suffi- 
ciently well consisted of a single elliptical and three circular 
Gaussian components. Each component (knot) is characterized 
by a flux density S at distance R from the center of the image 
with the relative R.A. and decl. equal to (0, 0), located at the 
flux density peak of the image, a position angle O with respect 
to the image center measured from north to east, and a size a 
corresponding to the FWHM of the Gaussian brightness 
distribution. Additionally, we performed modeling of the 
NRAO 530 ent- imaging results by fitting elliptical Gaus- 
sian components to the image-domain features. Table 2 gives 
the results of these three separate fitting procedures. To 
estimate uncertainties of components’ parameters, we have 
used an approach employed in high-frequency VLBI, e.g., 
Jorstad et al. (2005): the flux density uncertainties correspond 
to either 1096 of the flux density measurement or the noise level 
of the map obtained as the result of the modelfit task in 
DIFMAP depending on whatever is higher, while the position 
uncertainties are calculated as one-fifth of the synthetic 
CLEAN beam corresponding to the (u, v)-coverage for 
compact components, or equal to the beam size for diffuse 
knots with a size larger than the size of the beam. 

The source consists of two main components, CO and Cl, as 
indicated in Figure 5. The knots are separated by ~30 pas, with 
CO being a factor of ~3 brighter than C1. We associate CO with 
the core of NRAO530 at 1.3mm and determine the jet 
direction on microarcsecond-scales to be ~—28° east of north. 
All images shown in Figure 4 indicate nontrivial structure in 
CO, elongated along the jet axis. The models super-resolve the 
substructure of CO into components COa and COb. In the case 
of DIFMAP, the model yields the best fit to COb as a point 
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Figure 4. Images of NRAO 530 from the 2017 April EHT observations produced using DIFMAP, eht-imaging, SMILI, DMC, and Themis, shown within the 
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Figure 5. Fiducial EHT total intensity image of NRAO 530; the color scale 
corresponds to the brightness temperature; the contours represent the total 
intensity, with a peak value of 1.4 x 10!? K, starting at 10% of the peak and 
increasing in steps of 10%; the circle in the bottom-right corner shows the size 
of the nominal EHT beam with an FWHM diameter of 20 pas. The C2 contour 
is shown with a dashed line as it is below the 10% level of peak total intensity 
in the averaged image. 


Table 1 
Characterization of the Fits to NRAO 530 Data 
Quantity DIFMAP eht-imaging SMILI DMC Themis 
X ME 0.974 0.638 0.858 0.517 1.211 
Xep 3.812 0.899 2.346 0.977 4.210 
unos 1.513 1.137 2.721 0.943 2.844 


Flux (Jy) 1.030 0.958 0.899 1.585 1.120 


Note. For the x? calculation, the original thermal error budget was inflated by a 
nonclosing systematic error amounting to 2% of the complex visibilities. 


source. We interpret the brighter component C0a as the VLBI 
core. The weaker component COb is located ~20 uas southeast 
from the C0a core, along the same jet axis as the COa and Cl 
components. Taking into account a small viewing angle of the 
jet, it is very unlikely that COb is a counter-jet component. 
Instead, we interpret COb as the result of the opacity 
stratification in the core, with COb representing the more 
optically thick part of the jet. The Gaussian component fits 


Table 2 
Parameters of Total Intensity Components 


Knot S (mJy) R (pas) © (deg) a (uas) T, (10'° K) 
a) (2) (3) (4) (5) (6) 


DIFMAP (Gaussian modeling) 


C0a 278 + 29 1.8 +4 85.6 +71 16.3 + 4.2 23 
COb 61 X 15 21.9 + 4 149.6 + 13 0.0 

Cl 124+ 15 35.0225 26.9 +9 15.0 + 3.8 1.3 
C2 106+53 57.5+20 5.2 + 24 33 + 16 0.2 
BG 871 1.8 46.8 1000 


eht-imaging (Gaussian modeling) 


C0a 210 0.0 et 10.6 4.4 
COb 74 18.8 154.2 5.0 6.9 
Cl 108 28.5 —21.2 15.8 1.0 
C2 158 44.3 —23.8 42.8 x 13.6 0.6 
BG 1000 iss I sa ie 


eht-imaging (image domain, no blur) 


C0a 228 0.0 13.8 x 12.3 3.1 
COb 86 18.9 154.8 10.3 x 11.5 1.7 
CI 110 35.8 —28.0 13.0 x 12.7 1.6 
BG 978 E 3o. she ae 


Note. The columns are as follows: (1)—designation of knot; (2)—flux density 
in mJy; (3)—distance from the core in uas; (4)—position angle with respect to 
the core in degree; (5)—size in microarcseconds; (6)—observed brightness 
temperature in 10? K. 


show the presence of C2 component, and although uncertain- 
ties of the parameters are large, there is a good agreement 
between the results of Gaussian fits of DIFMAP and eht- 
imaging modeling. 

We have calculated the observed brightness temperature of 
the knots, 7p, in the same manner as in Issaoun et al. (2022; no 
cosmological redshift or Doppler corrections accounted for). 
According to Table 2, CO has a higher Tg than the CI 
component, which supports the interpretation that it is a VLBI 
core at 1.3 mm. The observed brightness temperature of the 
core is close to the equipartition value of T.g~5 x 10? K 
(Readhead 1994), which implies an intrinsic (plasma frame) 
brightness temperature T'g = Tg(1 + z)/0 & Teg, if we 
assume the same Doppler factor of 6~9 on pas scales as 
estimated on mas scales (Jorstad et al. 2017). A robust 
comparison with Teq requires calculation of 7g either for an 
optically thin knot or for the flux density corresponding to the 
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Figure 6. A comparison between the leakage calibration D-terms estimated by Themis from the NRAO 530 data set and D-terms assumed by other imaging 
methods, following the procedures described in Section 2. Empty circles correspond to the LO band, crosses correspond to the H I band, and error bars represent the 


lo width of the Themis posteriors. 


turnover frequency in the spectral energy distribution (SED) of 
a knot. However, the spectral index obtained for NRAO 530 by 
Goddi et al. (2021) is steep, œa ~ 0.8 (S;o v ^), while the 
source is dominated by the core at 221 GHz, so that even for 
CO we can assume that we measure Tp either near the turnover 
frequency or at the optically thin branch of the SED. The 
brightness temperature is somewhat higher at lower frequen- 
cies, where the core is optically thick. This will be discussed in 
a separate paper devoted to multiwavelength observations of 
NRAO 530 (M. M. Lisakov et al. 2023, in preparation). The 
results shown in Table 2 also contain a large Gaussian 
component, BG, with a flux density of -1Jy, which 
characterizes the unconstrained source emission on larger 
angular scales, overresolved on inter-site EHT baselines. 


4. Linear Polarization Images 


We have performed LP imaging of NRAO 530 using four 
out of five algorithms described in Section 3: DIFMAP 
(CLEAN) eht-imaging, DMC, and Themis. SMILI 
implementation does not currently have the capability to 
reconstruct polarimetric images. Given the large and poorly 
constrained systematic uncertainties in the polarimetric calibra- 
tion of the EHT data at times when ALMA is not participating 
in observations, only 2017 April 6 and 7 data were used, except 
for DMC, which allows for a very flexible calibration. of 
complex polarimetric gains simultaneous with the image 
reconstruction (Pesce 2021), and hence could reliably fit the 
entire data set. Similarly, except for DMC and Themis, we 
disregarded JCMT single polarization data for the LP analysis. 


4.1. Consistency Test for D-terms 


For all methods other than Themis, we either fixed all 
polarimetric leakage D-terms following the multisource fitting 
procedures described in EHTC et al. (2021) and Issaoun et al. 
(2022), or we only solved for the SPT D-terms. The latter 
resulted in the SPT D-terms being no larger than 1%, with 
uncertainties comparable to or larger than the values them- 
selves, which justifies setting SPT D-terms to zero. In the case 
of Themis, we performed full fitting of D-term coefficients 
simultaneously with the image reconstruction as an additional 
consistency test. In Figure 6 we provide a comparison between 
the assumed and estimated D-terms. The results show a 
generally high level of consistency, particularly when con- 
sidering the difficulty of quantifying the systematics present in 
a fit to a single sparsely sampled data set. JCMT D-terms are 
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constrained very poorly as a consequence of the station 
operating with a single polarization receiver in 2017. There is a 
good correlation between the assumed and calculated D-terms, 
corresponding to a Pearson coefficient of 0.86 for the vector of 
64 D-terms (ignoring the uncertainties) and median residual of 
2% on a single complex D-term. Altogether, the consistency 
test results support fixing D-terms for the model fitting and data 
analysis. 


4.2. Linear Polarization Images 


Figure 7 presents LP images of NRAO 530 in 2017 April, 
obtained using a selection of methods. The absolute magnitude 


of LP (y Q? + U? for LP Stokes image-domain components Q 
and U) is indicated with white contours and the length of the 
electric vector position angle (EVPA) ticks, while the EVPA 
tick colors represent fractional LP magnitude. The Faraday 
rotation toward NRAO 530 is entirely negligible at 230 GHz 
(Goddi et al. 2021); hence, observed EVPA can be directly 
interpreted as an intrinsic source property. The grayscale 
background image and black contours correspond to the total 
intensity. Similarly as in Section 3.6, in Figure 8 we also 
present the method-averaged image, where the background 
image and cyan contours represent the LP magnitude, while the 
white contours describe the total intensity image map. Similarly 
as in the case of the total intensity images, we adopted images 
shown in Figures 7-8 to a common field of view of 128 uas. 
The contribution of the emission outside of this region is 
consistent with the image noise. 

Figure 8 serves the purpose of identifying the total intensity 
and LP components. As for the total intensity images, the LP 
images show similar structure across all algorithms, although 
there are some differences in the polarization properties of the 
jet features. The most consistent among the methods is the 
polarized feature PO near the core, which has a slight shift to 
the northeast with respect to the total intensity peak CO in the 
DIFMAP, eht-imaging, and DMC images. In the DMC 
image, the polarized feature in the southeast part of the C1 knot 
is most prominent, compared to images produced by the other 
methods. The DIFMAP and DMC images show a polarized 
feature, P2, associated with the diffuse northeastern total 
intensity component C2, which has the highest degree of 
polarization. A hint of this feature can be seen in the eht- 
imaging image, while the Themis image possesses a 
prominent polarized feature P1 that is rather a combination of 
polarization that can be associated with knot C1 and the 
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Figure 7. LP images of NRAO 530 produced using DIFMAP, eht-imaging, DMC, and Themis. The total intensity is shown in grayscale with black contours 
indicating 10%, 25%, 50%, and 75% of the peak LP intensity. Black dotted contours indicate 25%, 50%, and 75% of the peak polarized intensity. The ticks show the 
orientation of the EVPA, their length indicates linearly polarized intensity magnitude, and their color indicates fractional LP. Cuts were made to omit all regions in the 
images where Stokes Z < 10% of the peak brightness and P < 10% of the peak polarized brightness. The images are all displayed with a field of view of 128 pas, and 


all images are blurred to an equivalent resolution of 20 pas. 
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Figure 8. Schematic of the total-intensity (CO, Cl, and C2) and LP (PO, Pl, and 
P2) components in the EHT fiducial image of NRAO 530; white contours show 
the total intensity levels as in Figure 5; color scale and cyan contours represent 
the polarized intensity of the method-averaged image. 
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northeastern component. Discrepancies seen in the Themis 
image may be related to the spline raster resolution limitations 
imposed by the MCMC algorithm. We expect that the method- 
averaged LP image, shown in Figure 8, should reveal the most 
reliable image-domain polarization features. 


4.3. Linear Polarization Component Modeling 


The method-averaged polarized intensity image presented in 
Figure 8 can be described with three polarized features 
designated as PO, Pl, and P2. Similarly as in Section 3.7, we 
have modeled Stokes Q and UM visibilities in DIFMAP using the 
task modelfit, in eht-imaging Gaussian component geo- 
metric modeling framework, and we also used image-domain 
feature extraction methods to characterize the eht-imaging 
results. Table 3 gives parameters of the polarized components, 
which are as follows: (1) polarized component; (2) distance 
between the center of the component and the total intensity 
peak; (3) position angle of the component with respect to the 
center; (4) FWHM size of the component; (5) Stokes O 
parameter; (6) Stokes U parameter; (7) fractional degree of 
polarization; and (8) component EVPA. Since Q and U data 
were modeled separately in DIFMAP, the distance, position 
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Table 3 
Parameters of Polarized Components 
Knot R (yas) O (deg) a (pas) Q (mJy) U (mJy) m (%) EVPA (deg) 
(1) (2) (3) (4) (5) (6) (7) (8) 
DIFMAP (Gaussian modeling) 
PO 5.7 £2.1 47 + 32 15.3 + 0.7 9.8 + 4.0 14.8 + 5.0 5.2 + 1.3 28+7 
Pl 26.6 + 1.8 1+9 16.8 + 1.1 3.0 + 4.0 —7.7+5.0 6.7 +2.5 —34+ 11 
P2 55.8 + 1.0 6+5 14.0+ —8.4 + 4.0 13.8 + 5.0 18.1 + 4.2 6l +12 
Map E e. e —6.1 + 4.0 24.7 + 5.0 1.7+0.6 52 + 10 
ALMA —8.6 36.8 2.35 + 0.03 51.6 + 0.4 
eht-imaging (Gaussian modeling) 
PO 0.0 EL 10.6 9.5 10.1 6.6 24.4 
Pl 28.5 —21 15.8 11.8 —8.8 13.6 —18.4 
P2 44.3 —24 42.8 x 13.6 —6.9 13.0 15.2 59.0 
Map E er 13.6 17:2 4.1 25.8 
eht-imaging (image domain, blur 15 pas) 
PO 7.3 68 11.5 7.6 0.5 8.0 2 
Pl 31.1 12 10.4 5.6 —0.5 22 -3 
P2 54.5 7 20.4 —9.8 25,2 58 56 
Map —11.2 15.5 2.7 63 


Note. (1) Values presented for the Map region are obtained by integrating modes of Stokes Z, Q, and U in the entire VLBI field of view; (2) Parameters presented for 
ALMA are results obtained by Goddi et al. (2021) on 2017 April 6. The columns are as follows: (1)—designation of knot; (2)—distance from the core in 
microarcseconds; (3)—position angle with respect to the core in degrees; (4)—size in microarcseconds; (5)—Stokes Q value in millijanskys; (6)—Stokes U value in 
millijanskys; (7)—degree of polarization in percentage; and (8)—PA of polarization in degrees. 


angle, and size of a component are calculated as averages over 
the O and U values, and standard deviations of the averages are 
given as uncertainties. The uncertainties in Q and U flux 
intensity are equal to the rms of the Q and U map noise. Based 
on the uncertainties in R, O, and a values given in Table 3, 
there is very good agreement in the position and size of 
polarized components between Q and M4 images, which 
supports the robustness of the DIFMAP modeling. We identify 
polarized components PO, P1, and P2 with the total intensity 
components CO, Cl, and C2, respectively, with parameters of 
the latter given in Table 2. The distances of the associated total 
and polarized intensity knots agree within 1.50 uncertainties, 
and their sizes are almost identical, except the P2/C2 pair in 
DIFMAP modeling where the polarized component P2 occupies 
only a fraction (less than a half) of knot C2. However, the 
position angles of PO and P1 show some differences relative to 
those of CO and C1, respectively, which cause the centers of the 
polarized components to be placed at the eastern edge of the 
corresponding total intensity components. The other EHT 
Sgr A* calibrator, J1924-2914, does not exhibit a similar effect 
(Issaoun et al. 2022), despite the fact that the two sources were 
calibrated in the same way and the same D-terms were 
employed to correct for the instrumental leakage. Hence, we 
conclude that these shifts are intrinsic to NRAO 530 rather than 
caused by instrumental effects. 

The LP image in Figure 8 shows polarized component P2, 
which we associate with total intensity knot C2, detected with 
Gaussian component models. However, due to its diffuse 
nature, the surface brightness of the knot is comparable with 
the noise level, so that C2 is indicated as a marginal feature in 
Figure 5, yet the knot has the highest degree of polarization in 
the jet. As mentioned above, in the case of the DIFMAP model 
polarized emission of C2, knot P2, occupies only a part of the 
total intensity feature. If we fix the size and position of a total 
intensity component, C2, equal to those of P2, and search for 
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the total flux density corresponding to the region of P2, 
this will yield a low value of —27 mJy, which implies the 
fractional polarization in the P2 region as high as ~60%. 
A similarly high degree of polarization is obtained in the 
case of eht-imaging image-domain analysis. Therefore, 
the P2 region should have an almost uniform magnetic 
field, which implies synchrotron emission fractional polarization 
of Mma = (a + 1)/(a + 5/3) x 100— 7596. (Pacholezyk 
1970), assuming a spectral index for P2 of a~ 1. One of the 
possible explanations why the knot is prominent in polarization 
but not in the total intensity is that this feature could be 
caused by shear, ordering the field, but not by shocks that 
would make it bright in total intensity. For the geometric 
modeling of Q and U in eht - imaging, we assumed that each 
total intensity Gaussian component corresponds to an LP 
component with a constant EVPA and fractional polarization, 
obtaining a lower fractional polarization of C2/P2, which is 
nevertheless significantly larger than that of the core comp- 
onent PO. 

Image-domain feature extraction with eht-imaging was 
based on masking disjointed image regions by LP flux density 
following the procedures used in Issaoun et al. (2022). We 
obtained three components that could be identified with PO, P1, 
and P2 from Figure 8. The results given in Table 3 are mostly 
in qualitative, if not quantitative, agreement between the 
methods, although the parameters agree within the 1.5c 
uncertainties given for values obtained by DIFMAP. The 
differences that are present can be attributed to the low 
polarized intensity in the jet relative to the polarized intensity 
noise and the systematic uncertainties specific to different 
imaging methods. 

Visual inspection of Figure 7 suggests that EVPAs in the 
P0/CO region lie between 20° and 35° in all polarization 
images, except the eht-imaging image, which is also 
consistent with the result of the Q and UM map modeling in 
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DIFMAP and eht-imaging Gaussian modeling. The ent- 
imaging image shows the EVPA to be closer to 0°, which 
can be explained by the higher level of polarization noise in the 
image. All images show polarized emission in the P1 region, 
with the EVPA close to the jet direction, while the EVPA of P2 
is perpendicular to the inner jet position angle as seen in the 
DIFMAP and DMC images in Figure 7 and determined by 
modeling. 

We have summed all CLEAN modes within the DIFMAP 
total intensity image in Stokes parameter Z (including the flux 
density of the BG component; see Table 2), and Stokes 
parameters Q and MU. This has resulted in an integrated total 
intensity flux of 1.5 Jy. Polarization parameters integrated over 
the submilliarcsecond EHT image field of view are listed in 
Table 3 under the component name “Map.” These parameters 
agree reasonably well with those obtained by Goddi et al. 
(2021) for NRAO 530 at 221 GHz using ALMA data on 2017 
April 6 and 7. There is very little variability of the flux density 
or polarization parameters between the two days, with flux 
densities reported by (Goddi et al. (2021) equal to 
1.61 3:0.16Jy (April 6) and 1.57+0.16Jy (April 7), and 
polarization parameters obtained on 2017 April 6 reproduced in 
Table 3 under the component name "ALMA." Including Stokes 
parameters Q and U of the P2 knot in the calculation of Map 
values is essential for a good agreement with the polarization 
measurements obtained by Goddi et al. (2021). 


4.4. Constraints on Circular Polarization 


The circular polarization of NRAO 530 at 230 GHz observed 
by ALMA (but unresolved) is consistent with zero (Goddi et al. 
2021). The high resolution of the EHT observation can be 
favorable for a detection of the circular polarization in the jet of 
NRAO 530, which is important for determining the plasma 
composition (Wardle et al. 1998). To find whether there is a 
statistically significant signature of spatially resolved circular 
polarization in the EHT VLBI data, we employed the same 
procedure as Issaoun et al. (2022), based on exploring the 
posterior distributions of the DMC fits to the observational data. 
Using 1000 images drawn from the posterior distribution, and 
evaluating the mean and standard deviation of circular 
polarization in each reconstructed pixel, we find no pixel with 
a Stokes V detection of significance larger than 1.4c. Hence, 
we conclude that we did not find statistically significant circular 
polarization resolved with the EHT resolution of ~20 pas. For 
comparison, the same procedure applied to the LP maps 
confirms detection at the ~7o level. Because of large image 
reconstruction uncertainties, this statement does not translate to 
a robust upper limit on the fractional circular polarization. 


5. Discussion 


The high-resolution total and polarized intensity images of 
the quasar NRAO 530 at 230 GHz obtained during the EHT 
campaign in 2017 April indicate a jet direction of ~—28° on 
pas scales. This jet direction is different from those seen at 
86 GHz (PA~ 10°; Issaoun et al. 2019) and 43 GHz 
(PA ~ —3°; Weaver et al. 2022), implying a wiggling of the 
jet with the distance from the core. In fact, the PA of 
component C2/P2 (see Table 3) agrees with that of the jet 
direction at 86 GHz within the uncertainties, which can be a 
sign of jet curvature already at this location. If that is indeed the 
case, the EVPA of the P2 component could be more aligned 
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with the local direction of the bending jet, which supports the 
origin of the emission in the shear. The curvature could be 
caused by an imbalance between the pressure inside and 
outside the jet resulting in the development of instabilities in 
the flow, or it could be connected with jet precession. An 
analysis of the jet structure from subparsec to kiloparsec scales 
using multiwavelength VLBI images of NRAO 530 contem- 
poraneous with the 230 GHz EHT image will be presented in 
Lisakov et al. (2023, in preparation). 

The brightest feature in the source, CO, which we associate 
with the VLBI core at 1.3 mm, is complex and consists of two 
components, COa and COb. Complex structure in the core was 
also revealed in the EHT images of the quasar 3C 279 (Kim 
et al. 2020). Variability of such core structures can be expected 
on short timescales with quasiperiodic oscillations produced by 
kink instabilities in the case of a magnetically dominated jet 
(Dong et al. 2020). The intrinsic brightness temperatures of 
both CO components are lower than the equipartition brightness 
temperature (see Table 2), even if the Doppler factor of the jet 
on microarcsecond scales is lower than that found on 
milliarcsecond scales. A possible explanation is that the energy 
density of the jet on microarcsecond scales is dominated by the 
magnetic field. 

Homan et al. (2006) analyzed the brightness temperatures of 
VLBI cores at 15 GHz of a large sample of AGNs. They found 
that, in a quiescent state, cores have a narrow range of intrinsic 
brightness temperatures close to the value expected for 
equipartition, while during outbursts the characteristic intrinsic 
brightness temperature significantly exceeds Teq. The authors 
have estimated that in active states, the energy in radiating 
particles exceeds the energy in the magnetic field by a factor of 
~10°. This implies that somewhere between the 230 and 
15 GHz cores, there is a region where a magnetically 
dominated jet changes its state and reaches equipartition 
conditions. Moreover, Jorstad et al. (2007) suggested, based on 
quasi-simultaneous multiepoch observations at 43 GHz 
(VLBA), 86 GHz (BIMA ^^, and 230/345 GHz (JCMT), that 
this region is located between the 86 and 43 GHz VLBI cores, 
while Lee (2013) found that 86 GHz VLBI cores in a sample of 
compact radio sources observed with Global Millimeter VLBI 
Array have a characteristic intrinsic brightness temperature 
lower than the equipartition temperature. The observed low 
brightness temperature of CO somewhat supports this scenario. 
However, analysis of the simultaneous VLBI observations at 
230 and 86 GHz is needed to draw firm conclusions. 

The polarized features in the jet are shifted with respect to 
the total intensity components with which they are associated. 
The shifts are almost transverse to the inner jet direction—the 
jet direction is from southeast to northwest, while the shifts are 
to the east. This feature appears robust in the image 
reconstructions. One viable interpretation of this phenomenon 
involves a difference in ordering of the magnetic field 
transverse to the jet, with the magnetic field on the eastern 
part of the jet being more ordered than that of the western part. 
This can be understood if the magnetic field in the jet consists 
of two components—ordered and turbulent, with the turbulent 
component being more pronounced on the western side of the 
jet, while a helical magnetic field can represent the ordered 
component. Indeed, a helical structure of the magnetic field is 
supported by the EVPA behavior in different polarized 
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components. The EVPA undergoes significant changes along 
~60 uas of the jet length: it goes from being oblique in the core 
to aligning with the jet direction in P1 ~30 uas to perpend- 
icular to the inner jet direction in P2 ~55 uas from the core. In 
addition, a higher polarization at the western side can be 
attributed to velocity shear between the jet boundary and the 
ambient medium being stronger on the western side, which is 
consistent with the explanation of the nature of component P2/ 
C2. However, the EVPA behavior is more aligned with a 
helical magnetic field, although we do not have a sufficient 
resolution to resolve the EVPA structure across the jet. Recent 
work by Pushkarev et al. (2022) analyzed the polarization 
properties in a large sample of AGNs observed with the VLBA 
at 15 GHz within the MOJAVE survey, which presented 
transverse fractional polarization slices beyond the core region 
for 307 sources. The authors have found that the majority of the 
sources (including NRAO 530) show a clear increase in 
fractional polarization toward the jet edges, with an asymmetric 
profile. They suggest that the observed patterns of polarization 
can be explained by helical magnetic fields with different pitch 
angles and different geometric and kinematic parameters of the 
jet, which is surrounded by a sheath. It therefore appears that a 
helical magnetic field structure persists in the quasar jet from 
subparsec to hundred-parsec scales. 

In addition to EHT data, we have obtained optical 
polarization observations of NRAO 530 several days before 
and after the EHT campaign. Figure 9 shows polarization 
measurements in the optical R band, at 221 GHz for the 
arcsecond resolution of ALMA (Goddi et al. 2021), and those 
of the EHT resolved components at 230 GHz. The optical 
observations were carried out with the 1.8 m Perkins telescope 
(Flagstaff, AZ) using the PRISM camera. The description of 
data reduction can be found in Jorstad et al. (2010). When 
comparing to the optical and ALMA values, the best agreement 
across wavelengths is found with the P2 component, indicating 
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dominance of the magnetic field direction perpendicular to the 
microarcsecond scale jet direction in the source. This suggests 
that the magnetic field direction in the optical emission region 
is similar to that of P2, which has the highest degree of 
polarization in the jet. 


6. Summary 


Using the 2017 EHT observations, we have obtained the first 
total and polarized intensity images of the quasar NRAO 530 at 
230 GHz, with an unprecedented angular resolution of 
~20 uas. We have not detected statistically significant 
variability of the source structure on a timescale of 3 days, 
and hence we have combined the data over epochs and 
frequencies to generate final images. We have used a variety of 
methods to construct the total and linearly polarized intensity 
images and averaged the final reconstructions produced by 
different methods to obtain the fiducial total and polarized 
intensity maps. 

On microarcsecond scales, the quasar has a bright core of 
7300 mJy, with a jet extending up to 60y:as along a 
PA ~ —28°. The core has substructure consisting of two 
features, COa and COb. We associate the 1.3mm VLBI core 
with the brightest knot, COa, while COb, located at a 
PA ~ 4-155? and —3 times weaker than COa can be a more 
opaque region of the core. The core is polarized, with a degree 
of polarization of ~5%—8%. The extended jet includes two 
features, CI/Pl and C2/P2, with mutually orthogonal 
directions of polarization, parallel and perpendicular to the 
jet, respectively, suggestive of a helical magnetic field 
structure. The centers of the polarized knots PO and P1 are 
shifted slightly transverse to the jet with respect to the centers 
of the corresponding total intensity knots CO and Cl. This can 
be caused either by a more ordered magnetic field on the 
eastern side of the jet, connected with a helical structure of the 
magnetic field, or a stronger interaction between the jet and the 
external medium on this side. 

Feature C2 has the highest degree of polarization in the jet, 
~20%. However, polarized knot P2, which occupies only a 
fraction of the total area of C2, has a peak polarization of 
~60%. The latter implies a nearly uniform magnetic field in the 
P2 region. In fact, the optical polarization measurement nearest 
in time to the EHT campaign has an EVPA similar to that in 
P2. NRAO 530 is a bright ^-ray source, and the origin of high- 
energy emission in blazars continues to be highly debated. We 
can speculate that the P2 region, with a very uniform magnetic 
field, might be a location where nonthermal optical and high- 
energy emission originates. Future multiepoch and multi- 
wavelength observations of the source simultaneous with EHT 
campaigns can provide a significant advance toward under- 
standing the location of the high-energy dissipation zone in 
blazar jets. 
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